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I. INTRODUCTION 

Quantum entanglement is one of the most remarkable 
features of quantum mechanics and it leads to powerful 
applications like quantum cryptography, dense coding 
and quantum computing [l|, 0] ■ 

It is well known that it is very hard to check whether 
a given density matrix describing a quantum state of 
the composite system is separable or entangled. There 
are several operational criteria which enable one to de- 
tect quantum entanglement (see e.g. [S] for the recent 
review). The most famous Peres- Horodecki criterion 
[3, 0| is based on the partial transposition: if a state 
p is separable then its partial transposition (]l(g)r)/5 is 
positive (such states are called PPT state). The struc- 
ture of this set is of primary importance in quantum in- 
formation theory. Unfortunately, this structure is still 
unknown, that is, one may easily check whether a given 
state is PPT but we do not know how to construct a 
general quantum state with PPT property. 

Recently Q we proposed a large class of bipartite 
PPT states which are based on certain cyclic decompo- 
sition of the total Hilbert space (8) C'' — we called 
them circulant states. The crucial property of this class 
is that a partial transposition of the circulant state 
has again a circular structure corresponding to another 
cyclic decomposition of C'^C'*. Interestingly, many 
well known examples of PPT states fit the class of cir- 
culant states 0. 

In the present paper we generalize the construction 
of circulant states to multipartite systems — N qudits 
living in (Xi . . . C {N copies) . This space can be 
in a natural way decomposed into subspaces of 

dimension d or d subspaces of dimension d^^^. Multi- 
partite circulant state is defined as a convex combina- 
tion of positive operators supported on these orthogo- 
nal subspaces. It turns out that the family of partial 



transpositions maps these states into circulant opera- 
tors supported on another family of subspaces related 
by a circular multipartite structure. Again, we show 
that many well known examples of multipartite states 
belong to our class. 

Recently, there is a considerable effort to explore 
multipartite systems [6|-(l3j and multipartite circulant 
states introduced in this paper may shed new light on 
the more general investigation of multipartite entangle- 
ment. 

The paper is organized as follows: for pedagogical 
reason we first illustrate our general method for d — 2 
in Section ini We recall basic construction for 2 qubits 
from 5] and then analyze in details 3 qubit circulant 
states. We illustrate our construction with well known 
examples of 3 qubit states from the literature. Then 
we present a general construction of N qubit states. 
Section Hill discusses circulant states of N qudits. Final 
conclusions are collected in the last section. 



II. AT-QUBIT STATES 

A. 2 qubits 

Consider a density matrix living in (g) which is 
given by 

P = Po + Pi , (1) 

where po and pi are supported on two orthogonal sub- 
spaces 

Eq — spanjeo (8i eo , ei (g) ei} , 

El = span{eo®ei ,ei(g)eo} , (2) 

and {eo,ei} is a computational base in C^. It is clear 
that {Eo,Si} defines the direct sum decomposition of 



C2(g)C2, that is 



where the matrices a = [oy] and b — [bij] read as follows 



So® Si = C^^C^ . 

We call it a circulant decomposition because its struc- 
ture is determined by the cyclic shift S' : ^ 
defined by 



Sci^ e^+i , (mod 2) 



One finds that 



(]l®5)Eo , 



(3) 



(4) 



and hence 



Po ^ 



Pi = 



i,J=0 



= (10^) ^ 6ye„ ®ey (l^Sy , (6) 

\hj=0 J 

where Cij := |ei)(ej|, and one adds mod 2. Now, since 
Pq and pi are supported on two orthogonal subspaces 
So and Si one has an obvious 

Proposition 1 p defined in (QJ) is a density matrix iff 

• a = [aij] and b = [bij] are 2x2 semi-positive 
matrices, and 



Tr(a 



1 



Now, the crucial observation is that partially trans- 
posed matrix p'^ = (11(g) r)/9 belongs to the same class 
as original p 



( 



P = 









aoi 












bw 


foil 




aio 






an , 



(7) 



and 



P = 



aoo 






aoi \ 




boo 


boi 






bio 


bii 




aio 






an / 



(8) 



aoo ^01 
bio ail 





boo 


aoi \ 




\ aio 


foil / 



(9) 



that is, both p and p'^ are circulant bipartite operators. 
Therefore, one arrives at 

Proposition 2 A circulant state represented by (QJ) is 
PPT iff a — [oij] and b — \bij\ are 2x2 semi-positive 
matrices. 

Note, that matrices a and b may be rewritten in the 
following transparent way 



(■5^ and similarly 



a ~ a ol -\- b o S , 



b = b ol-\- a o S , 



(10) 
(11) 



where x o y denotes the Hadamard product of two ma- 
trices X and y. 

Examples. 1. Bell states: \^^) = (|00) ± |11))/V2 



/ 1 • 



and for = (|01) ± |10))/\/2 



(12) 



1 / 



2. Werner state \\A 





/. . 




1 


• 1 


±1 • 


2 


• ±1 


1 • 




I- • 





(13) 





fl-p ■ 




1 


l+p 


-2p ■ 


4 


■ -2p 


l+p ■ 









(14) 



with "1/3 < p < 1- PPT condition implies well known 
result p < 1/3. 

3. Isotropic state [l^ 



1 = 



/l 



V 2p 



l-p 



2p \ 



l-p 



(15) 



2 



with — 1/3 < p < 1. Again PPT condition implies well 
known result p < 1/3. 



and for any two binaries fi and v define 



4. 0(2)® 


0(2)-invariant state 






= (I 




8)5'') A 








/ a + 2b 


2b- 


-a] 




One easily finds 










-J 


a + 2c 


a — 2c 




, (16) 


Aoi = span{eo<5 








eoj 


a -2c 


a + 2c 




5 eo <8 


ei , ei 


5 ei (8 




\2b-a 


a + 


2b ) 




Aio — span{eo$ 


5 61 8 


eo,ei(? 


5 eo (8 


ei} 


with a,b,c > and a + b + 


c = 1. It is 


clear that O is 


All — span{eo$ 


5 61 8 


ei , ei 


5 eo (8 


eo} 


positive and O is PPT iff 








It is clear that 












- 2 ' 


1 

^^2' 




(17) 


7^23 = Aoo © 


Aoi 


© AioS 


3 All 





which reproduces well known result |16l | . 

B. 3 qubits 

There are in principle two ways to generalize the 
2-qubit circulant decomposition for the case of three 
qubits. Either one decomposes 7^23 = 8) 8) 
into four 2-dimensional subspaces or into two 4- 
dimensional ones. 



1. 



20 20 2 02'- 



Let us define 2-dimensional subspace 

Aoo = span{eo (8 eo (8 eo , ei (8ei (8ei} (18) 
I 



(19) 



(20) 



(21) 



Now, we construct 3-qubit density matrix p of the fol- 
lowing form 

P = Poo + Poi + Pio + Pii , (22) 
where each is supported on A^^ . One has therefore 



Pm- = ""^r^ ® S^e^.S^'* (8 S'^e.jS"* , (23) 

which generalizes 2-qubit construction ^ . Positivity of 
p is guarantied by positivity of each 2x2 matrix [x'^'"''], 
and the normalization Trp = 1 is equivalent to 

Tr + + -f = 1 . (24) 

One obtains therefore the following block matrix 



■'-oo 

•^00 






•''01 

•''01 




(10) 
^00 

(11) 

^00 


^(10) 
•''01 

(11) 






(11) 

•''10 


(11) 

Xii 




(01) 

■^10 






(01) 
Xii 

• 4ry 



(25) 



where vertical and horizontal lines remind us about the 
splitting into blocks corresponding to tensor product 
structure (8 (8 . Double lines introduce split- 



ting into four blocks corresponding to (E) (C^ ® C^) 
and then single lines introduce splitting within each 
4x4 block corresponding to the second tensor product 
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(C^(8)C^). Now, let us perform partial transposition. 
There are three independent transformations 



Examples. 1. GHZ state [Ij 



with binary indices a and /3: 

Toi = t(g)t(g)T, 
Tio = ]1(X)T(8) 1 , 
Til = ]l(g)T(g)T . 



(26) 



(27) 



Now, it is easy to see that Ta^p has the same circular 
structure as original p defined in (|60p with new 2x2 
matrices that is, 

1 1 

i,j—0 

(28) 

Proposition 3 T/ie sei o/ 2 x 2 matrices is 
given by 

wif/i addition modulo 2. 

Using straightforward definition that p is (a/3)-PPT if 
''"q/3 P > 0, one has the following 

Theorem 1 A 3-quhit circulant state p is (a/3) -PPT iff 

for all binary p and v. 

It is clear that ([^5]) generalizes 2-qubit circulant state 
d?]). Note, however, that reducing 3-qubit state with 
respect to one subsystem one ends up with the following 
separable 2-qubit state 



Trip- t t 

fj,^i>—0 i,j—0 



(30) 



It is again circulant but has very special structure: [a^j] 
is diagonal with 

1 

a - V x^"'"^ 



and = 0. It is therefore clear that a general 2- 
qubit circulant state can not be obtained via reduction 
of (HSl). 



|GHZ) = — (1000) + 1111)) 



(31) 



does belong to circulant class which is easily seen from 
the corresponding density matrix 



1 1 
1 1 



On the other hand the well known W state 
|W) = ^(1001) + |010) + |100)) , 



(32) 



(33) 



is not circulant. The corresponding density matrix 
reads as follows 



/ 



pw 



V 



• 1 


1 • 


1 • 




• 1 


1 • 


1 • 




• 1 


1 • 


1 • 













(34) 



2. Bell states: the following 8 vectors 

= (-1)" (1(g) 5^ (g) S"')|GHZ) , (35) 

with binary q;,/3,7 define circulant states. These are 
3-qubit generalization of 2-qubit Bell states. Note that 

-0Q/37 e , (36) 

and the corresponding matrices x^^^^ read as follows 



^ \^ (-1)" 1 

3. Generalized 3-qubit isotropic state 



(37) 



1 - s 



®3 



s I GHZ) (GHZ I 



with s G [—1/7, 1]. One finds for x''^"' matrices 

^(00) ^ 1 / 1 + 3s 4s \ 

8 I 4s l + 3s / ' 



(38) 



(39) 
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and 



^(01) ^ ^(10) ^ ^(11) ^ L^i 



(40) 



The only nontrivial PPT condition comes from the pos- 
itivity of 



2-dimensional mutually orthogonal spaces A^^: either 
one introduces 



50 = Aoo ® All , 

51 = Aoi ® Aio , 



(43) 
(44) 



y 



mm 



1 M - s 4s 
8 I 4s 1 - ^ 



(41) 



which implies s < 1/5. Actually, it is well known [16 
that p is fully separable iff s < 1/5. 
4. 2-parameter 3-qubit state from [Toj : 



Aoo ' 
Aoi ' 



' Aio 
'An 



(45) 
(46) 



/I 



p(c, d) 



\1 



1 



• d 
c ■ 



d ■ 



1\ 



1/ 



(42) 



with —1/8 < c,d < 1/8. It was shown 19] that p{c,d) 
has positive partial transposes iff c = d. Moreover, this 
condition implies full separability. 



The construction is clear: 

50 = A^, , (47) 

fj.+i'=a 

51 = A,^, (48) 

whereas a second decomposition uses the following 
scheme 



A^,, 



(49) 



2. "8 = 404" 



There are two ways to construct 4-dimensional cir- 



Note, that using binary codes. So is constructed out of 
A^i^ with piy representing binary code for '0' (mod 2) 
and Si is constructed out of A^^ with fiiy representing 



culant decompositions of 3-qubit Hilbert space out of binary code for '1' (mod 2). One easily finds 



50 = span{eo ® Co ® Co , ei (g) ei ® ei , ei ® eo (8) eo , eo (8) ei (g) ei} 

51 = span{eo (8) eo (8) ei , ei (g) ei (g) eo , ei (g) Co (8) ei , eo (8) ei (g) eo} 



r 



(50) 



and 



10 = span{eo 8) eo 8) eo , ei (g) ei (g) ei , eo (8 ei (8 eo , ei 8) eo ® ei} 

11 = spanjeo ® eo ® ei , ei (8) ei (8) eo , eo (8) ei (8) ei , ei ® eo ® Co} 



(51) 
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It is clear that both decompositions are circulant, 
that is, 



and 



Si = (I(g)I«)S')I]o , 
Si = (I(g)I®S')So . 



(52) 
(53) 



Now we are ready to construct the corresponding 3- 
qubit circulant states: a circulant state corresponding 
to 



is defined by 



So © Si = , 



f = (To + fJl , 



with a a supported on Sq, . To define po and pi one needs 
4x4 matrices a and b. One has 



(54) 



1 1 
CTo = ^ ^ a^i 

i,j—Q 
1 1 

ci ^ ^ bf^i-^jefj_„®eij(g)SeijS* , (55) 

i,j—Q l-L^f—O 



where af^i-^^j and b^i.^j are matrix elements of a and 
b considered as matrices in the tensor product M4 = 
M2 (g) M2 = M2(M2). One obtains therefore the follow- 
ing block matrix 



aoo;00 

^00:00 


aOO:01 

&00;01 


000:10 

feQ0;10 


O00:ll 

&00;11 


^01:00 

aoi;00 


feoi;01 

aoi;Oi 


^01;10 

001:10 


&01;11 

Ooi;ll 


aiO;00 

^10;00 


aiO;01 

&10;01 


0X0:10 

&10;10 


OiO;ll 

felO;ll 


^11:00 

ail;00 


&11;01 

ail:01 


&11;10 

011:10 


6ll;ll 

Oii:ii 



(56) 



Let us note, that the density matrix defined by ((25)) is 
a special case of (|56|) where the matrices a and b are 
given by: 



(57) 



is defined by 



and 



b = 



■^00 


^(11) 

■^00 


^(11) 

•^01 


^(00) 
•''01 


■^10 


(11) 

^10 


(11) 

^11 


^(00) 

•'ll 


^(01) 
■^00 


^(10) 

-'■oo 


^(10) 
•^01 


^(01) 
•''01 


■^10 


(10) 
^10 


^(10) 

•^11 


^(01) 

•'11 



(58) 



Similarly, a 3-qubit circulant state corresponding to 
So ® Si = H23 , 



e = Co + 6 



with supported on S^. It is defined via two 4x4 
matrices c and d: 



1 1 

1 1 

6 = X! X! d.i^;jveij®e^u®SeijS* . (59) 



One finds 



6 



Coo ; 00 

f^OOiOO 


C00:01 

f^OOiOl 


CoO;10 

rfoO;10 


CoO:ll ^ 

rf00;ll 


C01;00 

f^OliOO 


C01:01 

f^OliOl 


C01;10 

c^oiao 


C01;ll 

C^Ol;!! 


(^lOiOO 

ClO;00 


f^lOiOl 

ClO;01 


c'loao 

ClO;10 


C^IO;!! 

CiO;ll 


rfll;00 

Cll;00 


f^lllOl 

Cll;01 


c'liiio 

Cll;10 


dll;ll 

Cll:ll / 



(60) 



It is clear that partially transposed T^/sa and t^^^ have 
exactly the same circulant structure as original a and 
^. One easily finds for the corresponding partial trans- 
positions: 

1 1 

i,j—0 fi,}y—0 
1 1 

i,j—0 fi,}y—0 

where the new 4x4 matrices c|^(^^- and are given 

by: 

c["0l = c , = d , (62) 

and 

cliol = (]l®r)c , = (]l®T)d , (63) 

where 

[it(g>T)c]ii_,-j^ = Ci^.-j^ , [(]l«'r)(i]j^y„ = di^.j^ . 
Moreover, 

c[""l O (I® I) + O (5®!) , 

where 

1 = 1 + 5 = ■ (65) 

The similar formulae one easily finds for a'"'^' and b^"^^ 
defined by 

1 1 

1 1 

+ E E ^Sfi ^/^^ ® ® ^ ^y-^* • (66) 



Theorem 2 ^ circulant 3-qubit state a is (a/3) -PPT 
are semi-positive matrices. Similarly, 
a circulant 3-qubit state ^ is (a/3) -PPT iff c^°'^'^ and 
d^"^^ are semi-positive matrices. 

Let us observe that a 3-qubit circulant state may be 
reduced to the 2-qubit circulant state. Consider for 
example a density operator ^ defined in (|59p . Note that 
reduction with respect to the second factor gives 
1 1 
Tr2^o = E E 

1 

= ^ (Tr2c)ij Cij (g) Cij , (67) 

1 1 

Tr26 = E d,i^;3ti ^ip ® ^ ^ij^* 
1 

= {^^2d\je,j®Se,jS* , (68) 

and hence 

TraC = TraCo + TraCi , (69) 

is a 2-qubit circulant state. It is no longer true for the 
remaining reductions with respect to the first and third 
factors. One obtains 
1 1 

Tri^ = ^ ^ [c + d\i^-it, e^^ ® en , (70) 

1— O /x,i^— 

and 

1 1 

Tra^ = ^ Y [c + d]ip;ii/eii(g)e^^ , (71) 

2— fl^U—O 

which are not circulant states. 
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C. TV qubits 

Consider now a general case of N qubits living in 
(C^)®^. Again, there are two natural ways to de- 
compose the corresponding Hilbert space 7i2N : either 
into 2^~^ two-dimensional subspaces or into two 2^~^~ 
dimensional subspaces. 



'2" = 2 © 2 ( 



®2" 



Let us introduce a circulant decomposition of 7^2" 
into 2^~^ two-dimensional subspaces. Now each integer 
from the set {0,1,..., 2^-^ - 1} may be represented by 
a string of — 1 binaries (/ii . . . (J^n-i)- Let us define 
2-dimensional subspace 

A0...0 = span {eo (Xi . . . (Xi eo , ei (g) . . . (g) ei} , (72) 
and for any string of binaries {fii . . . fiN-i) define 

A^,...^„_, = ... (g5^"-OAo...o . (73) 

Introducing convenient vector notation 
M = (/^i, • • ■ ,MJV-i) , 

one has 

A^ = (]l®5'^)Ao , (74) 

with 

S"' (E) . . . (gS*^"-! , 
and Aq = A0...0. One clearly has 

H2« - A^ , (75) 

where the sum runs over all binary (A^ — l)-vectors fi. 

Now, let us construct a circulant A^-qubit state p 
based on ([75]) : 



(76) 



where each is supported on A^. One has therefore 



1 

E 



(]1«)S''^)* , 

(77) 



where [a;*^'*^] are 2x2 semi-positive matrices. Normal- 
ization of p implies 



Tra;''"' = 1 



(78) 



Now, partial transpositions are labeled by a binary (A^— 
l)-vectors cr — {ai, . . . ,aN~i) 

= ]l(gT°'i(g ... (gr'^"-! . (79) 

Note, that each partial transposition Ta-p belongs to the 
same class of circulant states 



with 

>1 = (]1®5'^) 



(80) 



{1(g) sf^y 



where the new 2x2 matrices y('^)['^l are given by the 
following formula 



(81) 
the 

(82) 



A state p is cr-PPT iff Ta-p > and hence one has 

Theorem 3 A circulant state p is cr-PPT iff y'^'^^^'^^ 
are semi-positive for all fi. 

Examples. 1. Generalized GHZ state 

|GHZ) = -^(|0...0) + |l...l)) , (83) 

does belong to circulant class which is easily seen from 
the corresponding density matrix 

x'^"^ ^ S,,,,.! . (84) 

2. Generalized Bell states: the following 2^ vectors 

V-a. = (-1)"(I®^")|GHZ) , (85) 

with a = 0, 1 and binary (A^— l)-vector v define circu- 
lant states. These are A^-qubit generalization of 2-qubit 
Bell states. Note that 

i^a. e A, , (86) 
and the corresponding matrices x'^'^"' read as follows 



1 (-1)" 
(-1)" 1 



(87) 
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3. Generalized iV-qubit isotropic state 



2. " 2™ = 2 



P = 



1 - s 



s I GHZ) (GHZ I 



with s e [-1/(2^ - 1), 1]. One finds for a;*'') matrices 
ro)_ 1 I l + (2^-l)s \ 



2^-is l + (2^-l)s / ' 



and 



1 - s 



(90) 



for fi ^ 0. The only nontrivial PPT condition conies 
from the positivity of 



which implies 



l-s 2^-is 
l-s 



s < 



(91) 



)N-1 



1 



The above condition guaranties full A^-separability of p 

4. 2-parameter A^-qubit state from [20|: for —1/2^ < 
c,d< 1/2^ one defines a set of matrices x^*^^ 



2^ I 1 1 / ' 



(92) 



for fj, corresponding to binary representation of 
{0,l,...,2^-2-l}, 



^ 1 c 
2^ I c 1 / ' 



(93) 



for /X corresponding to binary representation of {2^ ^- 
l,...,2^-i-l}, and 



2~ I d 1 



(94) 



for fi corresponding to binary representation of {2^~^— 
1, . . . , 2^ — 1}. It gene ralizes a 3-qubit state defined in 
(|42)) . It was shown [20] that the above A^-qubit circulant 
state has positive partial transposes iff c = d. Moreover, 
this condition implies full separability. 



There are several ways to construct 2^"-'^- 
dimensional circulant decompositions of A'^-qubit 
Hilbert space out of 2-dimcnsional mutually orthogo- 
nal spaces A^. The following choice 



So = , 
\^i\=o 

5^1 = A^, 
ImI=i 



(95) 
(96) 



where 

I /^l = /ii + . . . + , 

gives rise to the circulant structure 

Si = (]l®^-i®5)I]o . (97) 
Another construction goes as follows 



(a\k) 

A^, a = 0,l, fc=l,2,...,7V-l, 

(98) 

where the sum ^'-"l'^-' runs over all fi with /i^. — a. 
Note, that 2(^1^) displays circulant structure defined 

by 



)Af-fe-l\ " 



S(o|fc) . (99) 



We shall consider only one scheme with fc = iV — 1 and 
to simplify notation let us define 

Ea ■.= E^a\N-i) , a = 0,1, (100) 
which satisfies 

Si = (]l®^-i®5)So . (101) 

This very choice has clear interpretation: to define Sq, 
we sum over all /x = {fii . . . /ijv_i) which represent bi- 
nary code for a (mod 2). 

Now, let us construct a circulant state 

with cTq supported on Eq,. It is clear that 

1 N-2 
0-0=^^ aat:(3j (^ Ca^f)^ (g) Ci-j (g> Bi-j , (102) 
cx,/3i,j=0 k=l 
1 N-2 

a, 13 i,j=0 k=l 
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where oc and /3 are binary (A'^ — 2)-vectors with coor- 
dinates ttfe and respectively, and [aai;i3j] and [bai-f^j] 
are 2^~^ x 2^^^ semi-positive matrices. 
Similarly, one constructs a circulant state 

with supported on It is clear that 

1 N-2 

a,/3 i,j=0 /s=l 
1 iV-2 

q:,/3 iJ—0 k=l 

with 2^~^ X 2-^"^ semi-positive matrices [c,c«;j/3] and 

Now, let us consider partially transposed iV-qubit cir- 
culant operators. The corresponding partial transposi- 
tions are labeled by binary {N — l)-vectors 

To- := l^r'^'O ... ^r"""-^ . (106) 

Note, that both r<,.cr and have exactly the same 
circulant structure as original a and ^. One easily finds 
for the corresponding partial transpositions: 

where ^IT^ are again supported on S^: 

1 N-2 



c«,/3 i,j=0 
1 



N-2 



fe=l 



with the new matrices [c'^.^^] and [rf-^j.j^] which are 
defined by the following formulae: 



40] 



c , dl"! = d , 



(109) 



and 



Theorem 4 A circulant N-qubit state ^ is <t-PPT iff 
c[°'l and are semi-positive matrices. Similarly, a 
circulant N-qubit state a is <T-PPT iff aW and ftW 

are semi-positive m,atrices. 

Let us observe that an A^-qubit circulant state may be 
easily reduced to the (A''— L)-qubit circulant state (with 

N — L > 2): let li, . . . ,Il denote L distinct integers 
from the set {2, 3, . . . , iV — 1}. Then the partial trace 



Tr/j. jj^^ defines {N — L)-qubit circulant state with new 
^N-L-i ^ 2^-^-1 matrices c' and d' defined by 



c' = Tri,...;^c, d' = Tri,...i^d . 



III. JV-QUDIT STATE 



(112) 



Consider now the most general case of N qudits liv- 
ing in (C'')®^. Again, there are two natural ways to 
decompose the corresponding Hilbert space Ti^w: ei- 
ther into d^~^ — d-dimensional subspaces, or into d — 
-dimensional subspaces. 



d®d{ 



(110) 

where 7 is binary {N — 2)-vector and we treat c and 
as a matrices living in the tensor product M^^~^. 
It is therefore clear that r-^ c denotes the corresponding 
partial transposition of c in the tensor product 
Moreover, 

= d[^01o(I®I®^-2) + c[^O]„(5^|0iV-2-) _ (111) 



Let us introduce a circulant decomposition of TC^n 
into d^~^ two-dimensional subspaces. Now each integer 
from the set {0, 1, ... , d'^~^ — 1} may be represented by 
a string of A'^ — 1 dinaries (/ii . . . j-iN-i), i-c each /i^ € 
{0, 1, . . . ,d — 1}. Let us define d-dimensional subspace 

A0...0 = span {eo (8) . . . ® eo . . . , ed_i <^ ... ^ea-i} , 

(113) 

and for any string of dinaries (/ii . . . /Ujv-i) define 

A^,...^„_, = (I®^''^® ... 5''"-OAo...o • (114) 
Introducing convenient vector notation 

H= (/Xi,...,/XAr_i) , 

one has 

= (]l®5'^)Ao , (115) 

with 

S^^ = S^'^^ ... ^S^"'-■' , 
and Ao = Ao...o- One clearly has 

H^N = A^ , (116) 
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where the sum runs over aU dinary {N — l)-vectors fi. 

Now, let us construct a circulant iV-qubit state p 
based on ([^5)1 : 

where each is supported on A^. One has therefore 



a;-7 e,: 



(118) 

where [x'^'^^] are d x d semi-positive matrices. Normal- 
ization of p implies 



J2 Tr-T^*") ^ 1 . 



(119) 



Now , let us look for the corresponding partial transpo- 
sitions Ta-p with Tct introduced in (I79|) . There is a crucial 
difference between qubit and qudit case: for qubits par- 
tially transposed state have exactly the same structure 
as the original one. It is no longer true for qudits. It 
was shown in [sl] that partial transposition gives rise to 
a new circulant structure governed by a certain permu- 
tation: let n be a d X d permutation matrix defined 

by 



for k = 1, . . . ,d—l. It turns out Q that partially trans- 
posed matrix Ta-p is related to the following circulant 
structure: 



Al^l = (]l®5'^)A[rl , 



where 



and 



A|r' = (]l®n'^)Ao , 



One clearly has 



7^,«=0 aw , 



(121) 



(122) 



(123) 



for each binary {N — l)-vector cr. 

One finds therefore the following cr-circulant struc- 
ture for Ta-p 



At 



(124) 



Ileo = eo , Ilefe = e^-fc 



(120) 



with 



{t®si^){t®\r) 



(125) 



where the new d x d matrices are given by the 

following formula 



d-l 



^y- ^[M+fe^]o(n-S''=) . (126) 



fc=0 

For d = 2 one finds 11 = I and the above sum reduces 
to two terms, only. One therefore recovers (|82p. 

Theorem 5 A circulant state p is cr -PPT iff 
are semi-positive for all dinary fi. 

Examples. 1. Generalized GHZ state 



^ d-l 

|GHZ) = — ^ Cfe® ... (^efc 



(127) 



k=0 



does belong to circulant class which is easily seen from 
the corresponding density matrix with 



where 



(128) 
(129) 



a=0 



11 



generalizes that is, lij — 1 for alH, j = 0, 1, . . . , d— 
1. 

2. Generalized Bell states: the following c?^ vectors 

i^a.. = (f^"'»5'')|GHZ) , (130) 

where the phase operator Q is defined via 

nck^aj'^ck, k = 0,l,...,d~l , (131) 

with uj — e^'^*/'^, define a circulant state for any a ~ 
0,1, . . . ,d ~ 1 and arbitrary dinary {N — l)-vector i/. 
These are iV-qudit generalization of A^-qubit Bell states 
Note that 



B. " = d^-i ffid^-i'' 



i^o.^ e , (132) 

and the corresponding matrices a:^*^^ read as follows 

x^f'^^S^^^h, (133) 
where the d x d matrix D, is defined by 

% = io^-' , (134) 

and generalizes a 2 x 2 matrix from (I87p . 
3. Generalized A^-qudit isotropic state 



P = 



1 - s 



s|GHZ)(GHZ| (135) 



with s e [—l/{d^ — 1), 1]. One finds for a;*^'^-' matrices 
^(0) _ ^ / 1 + (d^ - 1). d^^^s ] 



and 



M - 



l-s 
~d^ 



(137) 



for /X 7^ 0. The only nontrivial PPT condition comes 
from the positivity of 



l-s d^-\s 
d^-ifi l-s 



(138) 



which implies 



s < 



1 



d^-i + 1 



The above condition guaranties full A^-separability of p 



There are several ways to construct d^~^- 
dimensional circulant decompositions of A^-qudit 
Hilbert space out of d-dimensional mutually orthogo- 
nal spaces A^. The following choice 

Sa= A^ , a = 0,l,...,d-l (139) 

— a 



gives rise to the circulant structure 



E„ = (]l®^-i®5")So 



Another construction goes as follows 



(140) 



(141) 



for a = 0, 1, . . . , d - 1 and fc = 1, 2, . . . , TV - 1. In the 
above formula the sum ^'^"''^^ runs over all fi with 
Ilk = O-. Note, that Sj-^i^-) displays circulant structure 
defined by 



) 2(o|fc) . (142) 



We shall consider only one scheme with fc = — 1 and 
to simplify notation let us define 

S„ := S(a|Ar_i) , a = 0,\,...,d- \ , (143) 

which satisfies 

S„ = (]l®^-i®5")So . (144) 

This very choice has clear interpretation: to define Sq, 
we sum over all /i. = (/ii . . . /iAr_i) which represent di- 
nary code for a (mod d). 

Now, let us construct a circulant state 



0- = ^ CTq 



Q=0 

with ctq supported on Sq,. It is clear that 

d-l Ar-2 

= E E (8) ® e,, ® S"e,,S"* (145) 

a, 13 i j=0 k=l 

where a and (3 are dinary (A^ — 2)-vectors with co- 
ordinates a/c and Pk, respectively, and [a^-^^^] is a set 
of d^~^ X d^~^ semi-positive matrices. This set gener- 
alizes two matrices a = a^*^-* and b = a'^-* in the qubit 
case, i.e. d = 2. 
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Similarly, one constructs a circulant state 

d-l 



a=0 



with supported on ^q,. It is clear that 



= E E ^S./3 ® (8) ^'^A ® ^"e^-'^"* (146) 

a,f3i,j=0 k=l 

with d^"-'^ X d^"-'^ semi-positive matrices [c|^l^^]. 

Now, each partial transposition To- gives rise to the 
new circulant structure: either 

SM=0A|:1, a = 0,1,..., d-l (147) 



with the cyclic property 
or 

with the same property, that is. 



(148) 



(149) 



(150) 



One easily finds for the corresponding partial transpo- 
sitions: 



d-l 



c(=0 



where are supported on 

d-l 



(151) 



a,/3 ij— 



fc=i 



with new set of matrices [c^^'!^^'] which are defined by 
the following formulae: 



and 



(152) 
(153) 



where 7 is binary [N — 2)-vector and we treat c^"-* ma- 
trices living in the tensor product M^^~^. It is there- 
fore clear that r-y c'-"^ denotes the corresponding par- 
tial transposition of c*^"' in the tensor product n-i ^ 
Moreover, 



d-l 



^c("+'')hO]„(n5/5^I»^-2) .(154) 



Theorem 6 A circulant N-qudit state ^ is cr -PPT iff 
are semi-positive matrices for a^O, l,...,o?— 1. 



IV. CONCLUSIONS 

We have constructed a large class of PPT states which 
correspond to circular decompositions of C'' (g) ... (g) C 
into direct sums of d- and (i^~^-dimcnsional subspaces. 
This class generalizes bipartite circulant states intro- 
duced in [5|. It contains several known examples from 
the literature and produces a highly nontrivial family 
of new states. 

There are several open problems related to this new 
class: the basic question is how to detect entanglement 
within this class of multipartite states. One may ex- 
pect that there is special class of entanglement witnesses 
which are sensitive to entanglement encoded into circu- 
lar decompositions, that is, circulant Hermitian opera- 
tors W e Mf^ such that 



Tr{Wpi € 
for all product states pi 



®Pn)>0, (155) 

. . (g) pn, and 

Tt{W() < , (156) 

for some circulant state ^. It is interesting to explore 
the possibility of other decompositions leading to new 
classes of multipartite states. Let us note, that so called 
W state of 3 qubits ([34]) does not belong to our class. 
Another important family of states which does not fit 
circulant class was introduced in [T^ : these are A'^-qudit 
states satisfying 



U (g) . . . (g)U p ^ pU(, 



(157) 



for all unitaries U € U{d). For = 2 it reduces to 
the Werner state |l3j which belongs to bipartite circu- 
lant class. However, it is easy to check that for A^ > 3 
states satisfying (|157p are not circulant. One may ex- 
pect the existence of other characteristic decomposi- 
tions which are responsible for the structure of sym- 
metric states governed by (|157p . Anyway, multipartite 
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circulant states introduced in this paper may shed new 
Hght on the more general investigation of multipartite 
entanglement. 
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